Evidence of the in the K^K pn reaction 
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Abstract 

The —>■ K~^K~pn reaction has been studied to search for the evidence of the 6"*" by detecting 
K'^K~ pairs at forward angles. The Fermi-motion corrected nK'^ invariant mass distribution 
shows a narrow peak at 1.524 it 0.002 + 0.003 GeV/c^. The statistical significance of the peak 
calculated from a shape analysis is 5.1 a, and the differential cross-section for the jn if~6+ 
reaction is estimated to be 12 zt 2 nb/sr in the photon energy range from 2.0 GeV to 2.4 GeV in 
the LEPS angular range by assuming the isotropic production of the in the ^n center-of-mass 
system. The obtained results support the existence of the Q'^. 

PACS numbers: 12.39.Mk, 13.60.-r, 14.20.Jn, 14.80.-j 
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I. INTRODUCTION 



Since the LEPS collaboration reported the observation of a narrow baryon resonance-like 
structure in the nK~^ invariant mass spectrum produced in ^yn —>■ K^K^n reactions a 
considerable number of experiments have been carried out to check the existence of the exotic 
baryon, now called the 0"^. The 0^ is a genuine exotic baryon with the minimum quark 
configuration of uudds, for which a narrow decay width and a light mass were first predicted 
by Diakonov, Petrov, and Polyakov using a chiral quark soliton model Although the 
LEPS result seemed to be supported by several experiments which reported positive evidence 
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the experimental 



for the existence of the G"*" in various reactions 
situation soon became controversial. 

Many experiments at the high energy, especially collider experiments, found no positive 
evidence in t 
statistics 
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le pKs invariant mass distributions with a good mass resolution and high 
" " p. 
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16|,ll7|,ll8|. A typical upper limit for the inclusive production rate 



for the G"*" is less than 1% of that of the A(1520). The production mechanism of the G"*" 
might be very different from those of ordinary baryons if the G"*" exists jiol . 

The CLAS collaboration searched for the G"*" in the ■yp K^n reaction in the photon 

energy range from 1.6 GeV to 3.8 GeV with an integrated luminosity of 70 pb~^ 20|. The 
upper limit for the 7p K^Q^ reaction was determined to be 0.7 nb. The non-observation 
of the G"*" might be explained by a weak K*NQ^ coupling [21I]. On the other hand, if the 
K* coupling constant is small, the photo-production cross-section of the G^ from a proton 



22|. 



could be much smaller than that from a neutron 

The experiment which is most relevant to the current study was also carried out by the 
CLAS collaboration 23| . The search was done by detecting all charged particles in the final 



state in 'yd pK^K^n reactions with one order of magnitude higher statistics than the 
previous experiment ^]. The neutron momentum was reconstructed by using the missing 
momentum technique, and the G"*" was searched in the nK~^ invariant mass distribution. No 
narrow peak was observed, and the upper limit (95 % CL) for the elementary 'yn K~Q~^ 
reaction was obtained to be ~ 3nb by using a phenomenological model based on the A(1520) 
production to estimate the probability that the spectator proton is re-scattered and gains 
enough energy to be detected by the CLAS detector. 



Other dedicated experiments using vr [2J], K [25|, and proton 



beams have also shown 
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no evidence for the G"*" production. Although the theoretical interpretation of those null re- 
sults contains some uncertainties due to model dependences in the cross-section calculations, 
a strong K*NQ~^ coupling is unlikely if the G"*" exists. 

A model independent determination of the G^ width is possible by investigating the 
reverse reaction of the G+ decay: K^n G+. The DIANA collaboration observed evidence 
of the G"*" in the pKs invariant mass distribution from K^Xe K^pX reactions in a bubble 
chamber. The G"*" width was estimated to be 0.36 ±0.11 MeV/c^ from the production cross- 
section '2^- This result is barely consistent with the 90%-CL upper limit of 0.64 MeV/c^ 
from the Belle collaboration obtained by analyzing events from secondary kaon interactions 
in the material of the detector 28[|. 

To summarize the situation, if the G"*" exists, 1) its production is highly reaction depen- 
dent, 2) the coupling to K *N must be small, and 3) the decay width must be less than 1 
MeV/c^. Thus, it is desirable to study reactions which are sensitive to the KNQ^ coupling. 
The quasi-free reaction 'jn —>■ K~Q^ is one of such reactions since a 7 can couple to a 
i^'"^i^'~pair. 

In this paper we present a study of the photo-production of the G"*" from a neutron by 
closely comparing it with the photo-production of the A(1520) from a proton in a deuteron. 
Because the LEPS detector has a symmetric acceptance for positive and negative particles, 
a similar procedure can be applied to the both analyses. The validity of corrections and 
event selection criteria can also be cross-checked. 

The analysis is performed using the data collected with the LEPS detector in 2002-2003, 
where the statistics has been improved by a factor of 8 over the previous measurement jll. 



II. LEPS BEAMLINE AND DETECTOR 



A photon beam in the energy range from 1.5 GeV to 2.4 GeV is produced at SPring-8 
by Compton back-scattering of laser photons from 8 GeV electrons in the storage ring. The 
energy of a scattered photon is measured by tagging the electron which is associated with 
the Compton scattering event by event. The energy resolution for the tagged photon is 10 
MeV, and typical beam intensity with a 351-nm Ar laser is 10^ photons/sec. The photons 
are alternatively injected into liquid deuterium (LD2) or liquid hydrogen (LH2) targets in a 
16-cm long cell made of aluminum. The height of the interior of the cell is 60 mm, and the 
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width is 40 mm at the entrance and 100 mm at the exit. The windows of the cell are made 
of a Kapton polyimide film of 125 fim in thickness. 

The LEPS detector is a forward magnetic spectrometer which consists of a start counter 
(SC), a silicon vertex counter, an aerogel Cherenkov counter (AC), three drift chambers, 
a dipole magnet, and a wall of time-of-flight scintillation counters (TOF). The aperture of 
the 0.7-T dipole magnet is 55 cm high and 135 cm wide. The pole length is 60 cm. The 
angular coverage of the spectrometer is approximately ±20 and ±10 degrees in the horizontal 
and the vertical directions, respectively. The distance from the SC to the TOF is 4 m. A 
typical momentum resolution, Ap/p, for a charged particle with p = 1.0 GeV/c is 0.6 %, 
and the TOF resolution is 140 ps. The details of the detector and the quality of the particle 



identification are described elsewhere 



29|. 



The event trigger requires a coincidence of signals from the SC and TOF. A particle with 
p < 0.3 GeV/c cannot reach the TOF. Signals from the AC are used to veto events with 
pair-created e~^e~ or a pion with p > 0.6 GeV/c online. A typical trigger rate is 90 Hz for 
the LD2 runs and 60 Hz for the LH2 runs. 

The integrated numbers of photons in the energy range 1.5 — 2.4 GeV and 2.0 — 2.4 GeV 
were 3.93 X 10^2 and 2.15 x 10^^ fo^ the LD2 target runs, and the corresponding numbers were 
2.52 X 10^^ and 1.34 x 10^^ for the LH2 target runs. The total numbers of events collected 
with the LD2 target and the LH2 target were 4.5 x 10^ and 2.2 x 10^, respectively. Events 
in the energy range 2.0 — 2.4 GeV have been used for the current study. 



III. EVENT SELECTION 



We select events of the type ^yd K'^K~X, where X denotes particles which are not 
required to be identified by the LEPS detector. 

The momentum of a charged particle is reconstructed from the track information, and the 
velocity is obtained from the track length and the TOF information. The mass of the charged 
particle is calculated from the reconstructed momentum and velocity. The momentum de- 
pendent mass resolution ctjv/ for a kaon is calculated by using the measured momentum and 
TOF resolutions. A charged particle is identified as a kaon if the reconstructed mass is 
within 3.5crM of the nominal value. Events with a K^K^ pair are selected, and the vertex 
point of the two kaon tracks is required to be within the target volume. A hit position of 



5 



a track in the vertical direction at the TOF wall is reconstructed from the charge ratio and 
time difference of the signals from both ends of a TOF counter, and the horizontal position 
is obtained from the counter segmentation. The consistency between the reconstructed hit 
position and the extrapolated track at the TOF wall is checked to remove events with a 
decay-in-flight kaon. A total of 25820 and 8675 events passed all the selection cuts for the 
LD2 runs and LH2 runs, respectively. 

The invariant K^K~ mass {M{K^K~)) distribution for the LD2 runs is shown in 
Fig. [T](a). A narrow peak at 1.02 GeV/c^ is due to — > K^K'^ decays. Events with 
1.01 GeV/c^ < M{K+K~) < 1.03 GeV/c^ account for approximately 74% of the K+K' 
events. The p{'~f, K^K^) missing mass {MM{'~f,K^K^)) distribution for the LH2 runs is 
shown in Fig. [Il^b). These events are dominated by elastic processes. The missing mass 
resolution for a proton is seen to be 10 MeV/c^. Inelastic events with a high missing mass 
value oiMM{-f,K+K~) > 1.08 GeV/c^ are 3.5 % of the selected K+K~ events. Events due 
to misidentification of a pion pair, which result in a low missing mass value, are estimated to 
be less than 1 %. The MM(7, K^K~) distribution for the LD2 runs is shown in Fig. [H^c). 
The struck nucleon in the initial state has been assumed to be at rest. The peak near the 
nucleon mass is wide due to the Fermi motion of the nucleon. Some of the events in the 
lower tail region are due to coherent processes [s^, which are identified as a small peak at 
1.88 GeV/c^ in the dlj, K~^K~) missing mass {MMdlj, K^K~)) distribution as shown in 
Fig. [Id). 

IV. MINIMUM MOMENTUM SPECTATOR APPROXIMATION 

Because the momenta of the target nucleons are not measured, some approximation is 
necessary to obtain the invariant mass of pK~ or nK~^ pairs from 'yd — > K^K~pn reactions. 
The processes of interest are sequential processes of quasi-free productions of A(1520) or O"*" 
and their decays; jp —>■ i^"'"A(1520) —>■ K^K~p and 7^ K^Q^ K^K^n. We call the 
remaining nucleon which is not associated with the quasi-free processes a spectator. The 
spectator momentum due to the Fermi motion is approximately ~80 MeV/c for a deuteron. 
And it is small compared with momenta of a photon and kaons which are detected by the 
LEPS spectrometer. Therefore, the simplest approximation is to neglect the existence of a 
spectator. In this case the nK~^ invariant mass is obtained by calculating a (7, K~) missing 
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mass {MM {'J, K~)) with the assumption that the struck neutron is at rest in the initial state 
and always on-shell. We call this approximation the free nucleon approximation (FNA). A 
Monte-Carlo simulation study shows the mass resolution of the 0"*" using the FNA is ~30 
MeV/c^, which is mainly determined by the Fermi motion of a neutron. 

The minimum momentum spectator approximation (MMSA) has been developed in order 
to improve the mass resolution. In the MMSA a spectator is assumed to have the minimum 
momentum for the given total energy-momentum = (-Epn, ~Ppn)) of a pn pair, which 
is in turn assumed to be equal to the missing energy-momentum of the — > K^K~X 
reaction: 

Ppn = Pmiss ^P^+Pd- Pk+ - Pk- ■ (l) 

This assumption is not valid for inelastic events with an additional pion. However, these 
events can be removed easily as it will be shown below. Note that is derived from 
measured quantities and the deuteron mass {pd — (md,0)). The magnitude of the nucleon 
momentum (pcm) in the pn center-of-mass system is then given by 

•>/ (Mp„ mp rrin) (Mp„ - nip + rrin) {Mp^ + mp- rrin) {Mp^ -nip- rrin) 

in terms of the proton mass (nip) , the neutron mass (m„ ) , and the invariant mass of a pn 
pair (Mp„ = ). If, for a particular event, Mp„ is found to be less than nip + nin due to 
finite detector resolutions and coherent contributions, it is set equal to nip + m„. 

The momentum of a nucleon in the laboratory system has the minimum magnitude if the 
direction is anti-parallel to that of the total missing momentum. This topology is assumed 
in the MMSA. The minimum momentum, Pmin, is defined as the component of the spectator 
momentum in the direction of the missing momentum. Thus, it is given by 



Pmin = —pCM ■ 



1/2 I o \ P miss] /<^\ 

+ ypcM + mN- ^ — > (3) 



Mpn ^ Mpn 



where m^v is the mass of a nucleon which is assumed to be a spectator. With this approx- 
imation, the momentum component of the other nucleon in the direction of the missing 
momentum is given by 

Pres I P miss \ Pmin ■ (4) 
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If we assume that a spectator is a proton, the momentum of a neutron in the final state is 
given by 

— ^ P miss /_\ 
P n = Pres ' T=y i" • (,<3j 

I P miss I 

The invariant mass of the nK~^{M{nK~^)) is calculated by using the above and a mea- 
sured momentum. The resolution for the O"*" mass using the MMSA is 11 MeV/c^, 
which is an improvement over the FN A by a factor of 3. 

Events which are not associated with quasi-free processes can be identified from a large 
bmini value. Coherent processes which have a deuteron in the final state are characterized 
by a positive Pmin value which is approximately equal to a half of \'p'miss\- Inelastic reactions 
which create a pion in addition to a kaon pair cause pmm to have a large negative value. 
Re-scattering processes cause the pmin distribution to be dispersed. By requiring |pmin| to 
be small, these background events can be reduced. 

The Pmin distribution for the selected K^K~ events is shown in Fig. [2](a). The main 
contribution from quasi-free processes results in a peak near zero. The contribution from 
coherent processes is seen as a bump near 0.15 GeV/c, and the inelastic events concentrate 
in the region below -0.1 GeV/c. The projection of the spectator momentum onto the axis of 
'p^miss (pf) is well approximated by pmin as shown in Fig. EJ^b) for a Monte-Carlo simulation 
of non-resonant K^K~ events. 

The MM{'j,K+K^) and MMd{'j,K+K-) distributions for events with Ipminl < 0.1 
GeV/c are shown in Fig. [3l The inelastic and coherent contributions are successfully re- 
moved with Ipminl < 0.1 GeV/c. 

Let s be the square of the total center-of-mass energy of the nK^K^ system obtained 
with the MMSA. The effective photon energy E^^^Hs then defined by 

Note there is a one-to-one relation between E^-I'f and s. The -E'^-^-'^ becomes close to 
when the magnitude of the Fermi momentum is small. For the events with a small E^^-l" , 
all of M{K~^K~), M{nK^), and M{nK~) have a small value close to a threshold. Since 
we do not identify the nucleon in the final state, A(1520) events and events from protons 
and neutrons may contribute in the small M{nK^) region in this case. Events with a large 
E^-l"^ value are also problematic because they are dominated by coherent events and events 
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with particle misidentifications. Therefore, in addition to the condition |j9min| < 0.1 GeV/c, 
events are required to satisfy the condition 2.0 GeV < E!^^^ < 2.5 GeV for further analysis. 
The value for E^f^ = 2.0 GeV is 2.15 GeV. Thus, the maxmimum Mnk is 1.65 GeV/c^ 
at the cut boundary. Events with E^-^-l" > 2.5 GeV have large Fermi momentum for which the 
MMSA is not a good approximation. The number of events with E^-^-l" > 2.5 GeV is small 
(658 events) compared to 14928 events with 2.0 GeV < E^-l'-f < 2.5 GeV. For the LH2 runs, 
we require events to satisfy 2.0 GeV < E^ <2.4 GeV and 0.9 GeV/c^ < MM(7, K+K~) < 
0.98 GeV/c^. A total of 6306 events have passed the requirements. 

In principle, m„ should be replaced by rrip in the case of the pK'^K" system. However, 
the difference between E!^^^ values calculated with m„ and rrip is less than 1 MeV. Therefore, 
we use the mean of rup and for the calculation of E^^^ . 

The Pmin distributions for events with 2.0 GeV < E!^^^ < 2.5 GeV are shown in Fig. HI 
Both the coherent and inelastic events are strongly suppressed, and the main peak due 
to quasi-free processes is well reproduced by a Monte-Carlo simulation for non-resonant 
'-fn —>■ K^K~n reactions using a realistic deuteron wave function [31]. The non-uniform 
structure in the higher tail of the distribution for the Monte-Carlo events is caused by the 
special treatment of setting Mp^ = nip + m„ when Mp„ becomes smaller than rrip + m„ due 
to the finite resolutions. 

In the Monte-Carlo study, the mass of a struck nucleon has been set to be off-shell so 
that the total energy of two nucleons in the center-of-mass system is equal to m^. The mass 
of a spectator nucleon is always set to be on-shell. 



V. RANDOMIZED MINIMUM MOMENTUM METHOD 

In this section, we develop a method to estimate the reasonable M{nK^) spectrum 
shape for background contributions by using only measured E-y and Ji^K- values. There is 
a strong correlation between pmin and MM(7, K~) for the signal Monte-Carlo events, while 
the correlation is very weak for non-resonant background events as shown in Fig. [5l The 
nature of background events is characterized by the absence of this correlation. Because 
Pmin of a background event has a random nature due to the Fermi motion, a reconstructed 
Pmin can be replaced by a computer-generated one without changing the shape of a M{nK~^) 
distribution for background events. 
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The first step of tfie randomized minimum momentum metliod (RMM) is to approximate 
tfie mass correction (AM = M{nK^) — MM{'-f, K^)) by a 2nd order polynomial function 
of Pmin (= AM'(pinin)) as showu in Fig. El^a). The quality of this approximation is quite 
accurate, and the standard deviation of AM — AM' is 4-5 MeV/c^ in the whole mass range. 
Fig. MJo) shows the M{nK~^) and MM{'y, K~) + AM' distributions for the signal Monte- 
Carlo events. It demonstrates that the mass correction is predominantly determined by Pmin, 
and other effects such as the directions of 'p^K+ and are small. The practical advantage 
of this simplified mass calculation is separation of input arguments for the mass function 
into two types: one which depends on only E-y and "p^x- , and the other which also depends 
on 'j) K+- The original M{nK^) with the MMSA is a complicated function of E^, ~P k-i and 
'j) K+- In the RMM, it is approximated by a function of MM{'y, K~) and Pmin- 

In the next step, the most probable M{nK^) spectrum shape for a given MMIj, K~) 
distribution is estimated by combining each MM{'y, value with randomized pmin val- 
ues for many times (10^ times in this analysis). In the generation, the Pmin distribution 
is assumed to have a Gaussian shape. Because there is a weak correlation between pmin 
and MM{'y, K~) near the tails of the MM{'y, K~) distribution, the mean of the Gaussian 
distribution must be varied as a function of MMIj, K^). This correlation is mainly caused 
by the difference between the kinematic domains of MM{'y, K~) and M{nK~^); the former 
can have a value below m„ + mK+, but the latter cannot. For the same reason, the standard 
deviation (a) of the pmin distribution must be varied near the mass threshold. The magni- 
tudes of changes in the mean and a are small compared to a typical a value of ~40 MeV/c^. 
Fig. M^c) shows the mean and ±lcr curves as functions of MM(7, K'). The same functions 
are used for the estimation of the M{pK~) spectrum shape from MM(7, K~^). 

The RMM is similar to the mixed event technique which is widely used for estimations 
of combinatoric background. Both methods require the independence of uncorrelated vari- 
ables, and a common problem is signal contamination. In the case of the RMM, the signal 
contamination causes enhancement of the background level in the region of interest. 

The final step of the RMM is to divide the real data events which are used for the seeds 
of the event generation into several sub sets. We call them seed sets. In the current analysis, 
events are sorted according to the value of M{nK'^) or M{pK~) which is most sensitive to 
the signal-to-noise ratio. The boundary for one of the seed sets is chosen to cover a signal 
region or more precisely a possible signal region. The events of the seed set for the signal 
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region would contain both signal events and background events. A good feature of the RMM 
is that the shape of the mass distribution generated with the signal events is very close to 
the shape of the mass distribution generated with the background events in the same seed 
set. Fig. Wi^) shows the RMM spectra generated with the G''" MC events and with the non- 
resonant K^K^ MC events both in the region of 1.50 GeV/c^ < M{nK+) < 1.55 GeV/c^. 
The difference in the spectrum shape is small. 

In the shape analysis in the following sections, the background spectrum is represented 
as a sum of several RMM spectra. If one of the seed sets is contaminated by signal events, 
the effect can be absorbed by reduction of the weight parameter in the summation of the 
RMM spectra. By using several RMM spectra for the shape analysis, some of the global 
inconsistencies caused by fluctuations in the seed distributions and incomplete treatment of 
the correlation between pmin and MM(7, K^) can also be compensated by small changes of 
the weight parameters. 

The significance of a signal contribution is calculated from the difference in log likelihood 
between fits with and without the signal contribution represented by a Gaussian function. 
Since the width is fixed to the value estimated by a Monte-Carlo simulation, the change in 
the number of the degrees of freedom is 2, which is taken into account for the significance 
calculation. 

The M{nK~^) distribution for (j) events selected with requirements of 1.01 GeV/c^ < 
M{K+K-) <1.03 GeV/c^ and 2.0 GeV < E^^^ < 2.5 GeV is shown in Fig. EKb). A fit 
to a spectrum generated with the RMM using all selected events with an equal weight is 
indicated by the dashed line. The solid line shows the fit results with three RMM spectra, for 
which the selected events are divided into the three seed sets according to the conditions: 
M{nK+) <1.5 GeV/c^, 1.5 GeV/c^ < M{nK+) < 1.6 GeV/c^, or M{nK+) > 1.6 GeV/c^. 
The log likelihood (-21nL) for the single RMM spectrum fit is 65.5 for the number of the 
degrees of freedom (ndf) equal to 51. It becomes to 50.4 for the fit with three RMM spectra 
for Andf = 2. Since the spread of a RMM spectrum from the seed set of 1 MeV/c^ width 
is larger than 30 MeV/c^, further increasing the segmentation for the seed sets does not 
improve the -21nL value more than Andf. 

The M{nK~^) distribution for the sum of 3000 non-resonant K^K~ MC events and 
300 MC events is fitted to a mass distribution consisting of three RMM distributions 
with seed regions of (I) M{nK+) < 1.50 GeV/c^ , (II) 1.50 GeV/c^ < M{nK+) < 1.55 
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GeV/c^, and (III) M{nK+) > 1.55 GeV/c^. The best fit, which is indicated by a sohd 
curve in Fig. [7](c), is obtained with the weight parameters of 0.651, 1.245, and 0.949 for 
the contrbutions from region (I), (II), and (HI), respectively. The -21nL value for the fit 
is 114.6 for ndf = 61. The -21nL value is improved to 58.4 for ndf = 59 by including a 
Gaussian function with a fixed width of 11 MeV/c^ to represent the contribution as 
shown in Fig. El^d). The statistical significance of the peak is calculated to be 7.2 a for 
A {-2lnL) / Andf = 56.2/2. The weight parameters become 1.14, 0.648, and 0.993 for the 
contrbutions from region (I), (II), and (HI), respectively. The sum of the RMM spectra with 
the fitted weight parameters, which is indicated by a dotted curve in Fig. Wlyd), reproduces 
well the original mass distribution (open circles) for the non-resonant K^K~ MC events. 
The signal yield estimated from the fit is 279 ± 36 events, which is consistent with the 
number of 0"^ MC events of 300. 

VI. ANALYSIS OF THE -fp /s:+A(1520) REACTION 

The dominant contribution in the selected K^K~ events is due to decays, which are 
rejected by the combined requirements of M{K^ K^) > 1.03 GeV/c^ and M{K^K^) > 
1.02 + 0.09 X {E^-^-l' - 2.0). The cut boundary is shown as solid lines in Fig. [8]^a). The energy 
dependent cut condition makes the signal acceptance more uniform than the constant cut 
condition as shown in Fig. Mjo). The M{nK~^) distributions for non-resonant K^K~ MC 
events before and after the exclusion cut are shown in Fig. [8]^c). The exclusion cut 
distorts the mass spectrum because the acceptance is high near the threshold where the 
momenta of the and are highly asymmetric. However, the mass dependence of the 
acceptance is not strong, and consequently the cut does not create a narrow peak. Note 
that the mass spectrum near the threshold is not affected by the exclusion cut. This is 
because the momenta of two kaons for events in the threshold region are highly asymmetric, 
which results in a high K^K~ mass. 

A total of 2078 events passed the exclusion cut, and the MMSA is applied to ob- 
tain M{pK~) by assuming the spectator is a neutron. The Dalitz plots {M'^{pK~) vs. 
M'^{K^K^)) before and after the exclusion cut are shown in Fig. [9l The A(1520) yield at 
M'^{pK^) ~ 2.3 GeV/c^ is higher in the lower M'^{K^K') region due to the LEPS detector 
acceptance. However, the events are not concentrated near the cut boundary. 
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Fig. [TOT a) shows the M{pK~) distribution. For the shape analysis, three RMM spectra 
are generated by setting the seed boundaries at 1.48 GeV/c^ and 1.56 GeV/c^ in M{pK~). A 
fit to the RMM spectra gives a -21nL value of 110.2 for ndf=58. The -21nL value is improved 
to 55.1 by including a Gaussian function with a fixed width of 16 MeV/c^ as the A(1520) 
contribution. The A(-21nL) of 55.1 for Andf=2 corresponds to a statistical significance of 
7.1a. The signal yield is determined to be 289±38 events from the fit. The fit result with 
the A(1520) contribution is represented by the solid curve. The dotted line is the sum of 
the RMM spectra, which represents the background. The dashed line shows a fitting result 
without the A(1520) contribution. 

The data points of the MM{'y, K~^) and MM{'y, K~) distributions for the LH2 runs 
are shown in Fig. [TUT b) as closed circles and open circles, respectively. The A(1520) peak 
becomes narrow because of no Fermi motion effect. No significant peak structure is observed 
in the MM{'y, K~) distribution. The mass distribution in the region of 1.47 GeV/c^ < 
MM{'y,K^) < 1.65 GeV/c^ is fitted to a second-order polynomial. The result gives - 
21nL=32.3 for ndf =27. An excess of the MM{'y,K+) yield over the MM{'y,K-) yield is 
seen near the NK mass threshold. The excess of the MM{'~f, K~) yield in the high mass 
regions is due to the reflections of the A(1520) events. The MM(7, K'^) distribution is fitted 
to a Gaussian function plus a second-order polynomial, and the A(1520) yield is estimated 
to be 143±17 from the fit. The LD2/LH2 ratio of the A(1520) yield is 2.02±0.36. This is 
consistent with the estimated ratio of 1.93 from the integrated numbers of incident photons 
and the target proton densities. 

VII. ANALYSIS OF THE -/n K 6+ REACTION 

For the analysis of the 7^ K^Q^ reaction, events with a 3rd charged track in addition 
to and K~ tracks are removed. This condition changes the total number of events from 
2078 to 1967. Most of the removed events are due to 7p —>■ K^K~p reactions with a neutron 
as a spectator, for which the LEPS detector has a finite acceptance to detect all the charged 
particles in the final state. 

Dalitz plots of M^{nK^) vs. M^{K^K~) before and after the (p exclusion cut are shown in 
Fig.[TTl No concentration of events near the cut boundary is seen. The M{nK^) distribution 
for the final candidate events is shown in Fig. [T2]^a). There is a narrow peak structure near 
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1.52-1.53 GeV/c^. The distribution is fitted to a mass distribution consisting of tliree RMM 
distributions witfi seed regions of M{nK^) < 1.50 GeV/c^ , 1.50 GeV/c^ < M{nK^) < 1.55 
GeV/c2, and M{nK+) > 1.55 GeV/c^ The -21nL value of the fit changes from 104.7 (for 
ndf =66) to 73.64 (for ndf=6A) by including a Gaussian function with the estimated width 
of 11 MeV/c^ to represent the signals. The statistical significance of the signal estimated 
from A(-21nL) is 5.2a. The peak position for the best fit is 1.524 ± 0.002 + 0.003 GeV/c^ 
where the systematic shift of the peak position by +3 MeV/c^ due to the MMSA and the </> 
exclusion cut is given as a systematic uncertainty. The signal yield is estimated to be 116 
±21 events from the fitted peak height and its uncertainty. The detector acceptance has 
been calculated by assuming the isotropic production of the B"*" in the 7n center-of-mass 
system, and the differential cross-section for the jn —>■ K^Q^ reaction is estimated to be 
12 ± 2 nb/sr in the LEPS angular range. 

There is a dip near 1.56 GeV/c^ even with the B^ contribution. However, with the current 
limited statistics, it is not clear if the dip is due to fluctuations or due to some interference 
effects. Since we assume the branching ratio B(B+ —>■ K^n)=0.5 in the calculation of the 
differential cross-section, possible interference effects between the signal and background 
amplitudes could result in a change of the estimated value. 

A fit of the M{nK~^) distribution to the mass distribution using a Gaussian function with 
a free width parameter has been carried out, and the best fit is obtained with a width of 
12.7 ±2.8 MeV/c^, which is consistent with the estimated width of 11 MeV/c^. 

Fig. fT2l(b) shows the Dalitz plot of M'^{nK~^) vs. M'^{pK~). Note a proton is assumed 
to be a spectator for the calculation of M{nK^) and a neutron is assumed to be a spec- 
tator for the calculation of M{pK~). The relatively large -21nL values for the fits of the 
M{nK~^) distribution could be due to the reflections of the A(1520) events which might be 
responsible for the remaining structure near 1.65 GeV/c^. To avoid a possible effect due 
to the reflection, we require events to satisfy M{pK^) > 1.55 GeV/c^ and restrict the fit 
region up to 1.65 GeV/c^. The fit qualities are improved, giving -21nL/n(i/=55.2/33 and 
-21nL/n(i/=24.8/31 for the cases with and without the B"*" contribution, respectively. The 
significance is unchanged because the change in A(-21nL) is small. Fig. [T3](a) shows the 
M{nK^) distribution after the A(1520) exclusion cut. 

To study the model dependence, we have varied the boundaries of the seed regions for 
the RMM spectrum generation: the narrow signal region case with the boundaries at 1.51 
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GeV/c^ and 1.54 GeV/c^, and the wide signal region case with the boundaries at 1.48 GeV/c^ 
and 1.57 GeV/c^. The fit results are essentially unchanged, giving a statistical significance 
of 5.2 (J for the narrow signal region case, and a significance of 5.1 a for the wide signal 
region case. Although the shape and magnitude of each RMM spectrum vary case by case, 
the resultant summed background spectrum for the best fit is similar to each other as shown 
in Fig. [TST b). Fine structures in the original M{nK~^) distribution compared to a typical 
mass spread of ~30 MeV/c^ due to Fermi motion are smoothed by using the RMM. The 
weak dependence of the fit results on the seed boundary condition is a consequence of the 
smooth nature of the uncorrelated background. The maximum difference in the fitted peak 
height with the various background models is approximately 5 %, which is much smaller 
than the fitting uncertainty of 18 %, and, therefore, neglected. 

We have also examined a fit to the mass distribution using a second-order polynomial to 
represent the background. Note the number of the fitting parameters for the fit is the same as 
the fit with three RMM background spectra. Fig. [H] shows the comparison of the fit results 
for the cases with the fitting regions of 1.43 GeV/c^ < M{nK^) < 1.65 GeV/c^ and 1.47 
GeV/c^ < M{nK+) <1.65 GeV/c^ Fit quality is always better with the RMM than with 
the polynomial background function. By using the polynomial background function, the 
-2lnL/ndf values for the wide fitting region are obtained to be 65.1/33 and 28.1/31 without 
and with the 0^ contribution, respectively. The A(-21nL) of 37.0 corresponds to a 5.7a 
significance. For the narrow fitting region case, the -2\nL/ndf values for the fits without 
and with the G"*" contribution are 58.4/27 and 23.1/25, resulting in a 5.6 a significance. 
For the same fitting region, the -2lnL/ndf values using the RMM are 51.2/27 and 21.2/25, 
giving a significance of 5.1 a. Thus, the statistical significances estimated from the fit results 
with RMM are smaller than those with the polynomial functions. The difference is caused 
by poor modeling of a background shape with the polynomial functions, especially without 
the 0+ contribution. The fit results with the various background models are summarized 
in Table [H The smallest significance of 5.1 a is considered as the 6"*" significance with the 
systematics taken into account. 

The validity of the statistical significance estimated from the A(-21nL) value is checked 
with 2 X 10^ sample mass distributions generated by a toy Monte-Carlo simulation program 
by assuming a spectrum shape for the non-resonant K^K~ events. The generated distri- 
butions are fitted to the mass distribution which uses the polynomial background function. 
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TABLE I: Summary of fit results that are used to study the significance of the @^ contribution. 



Background model 


A(1520) cont. 


Fit region 


-2lnL/ndf 


-2lnL/ndf Sij 


j;nificance 






(GeV/c2) without G+ 


with 0+ 




F?lVTA/r sDPrtrri rlpffinlt, sppd spts 


not, pyrin rlpfl 




1 04 7 /6fi 


73 64/64 


5.2 a 


RMM spectra, default seed sets. 


excluded 


[1.43,1.65] 


55.2/33 


24.8/31 


5.2 a 


RMM spectra, wide signal region. 


excluded 


[1.43,1.65] 


54.5/33 


24.3/31 


5.1 a 


RMM spectra, narrow signal region. 


excluded 


[1.43,1.65] 


55.9/33 


24.8/31 


5.2 a 


RMM spectra, default seed sets. 


excluded 


[1.47,1.65] 


51.2/27 


21.2/25 


5.1 a 


2nd-order polynomial. 


excluded 


[1.43,1.65] 


65.1/33 


28.1/31 


5.7 a 


2nd-order polynomial. 


excluded 


[1.47,1.65] 


58.4/27 


23.1/25 


5.6 a 



The difference of the -21nL values with and without a Gaussian function with the width of 
11 MeV/c^ is checked. The numbers of the samples with a significance of more than 4cr and 
5cr are 10 and 2, respectively. These numbers are consistent with the expected number of 
occurrences of high-significance samples due to statistical fluctuations. 

A background spectrum for MM(7, K^) and MM{'y, K'^) distributions can be obtained 
by using the RMM in a reversed way, where a measured M{nK^) {M{pK~)) and a ran- 
domized Pmin are used to simulate MM{'y, K~) {MM{'y, K~^)). The missing mass distribu- 
tions with a fit to the reversed RMM functions are shown in Fig. [151 Both MM(7, K~^) 
and MM{'y, K~) distributions are well reproduced by the background functions with - 
2lnL / ndf =79 .7 / 67 and 60.2/66, respectively. Since the -21nL values are not reduced by 
more than Andf by including a Gaussian function with a fixed width of 30 MeV/c^, the 
corresponding significances are less than 1. This demonstrates the importance of the nar- 
rowness of the width and the consistency between the measured and estimated values for 
the shape analysis. 

A photon energy independent (j) exclusion cut with the condition of M{K^K^) >1.05 
GeV/c^ is also tried, and the resultant M{nK^) and M{pK~) distributions are shown in 
Fig. [TBI Both distributions are well fitted to the mass distributions with the RMM and 
Gaussian functions with -2lnL/ndf ratios of less than 1. The peak positions have not been 
changed, and the signal yields for the G"*" and A(1520) are reduced by a factor of 25% and 
35%, respectively, compared to those obtained by the original exclusion cut. 
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Events with a 3rd charged track are examined to check if the narrow peak in the M{nK~^) 
distribution is due to quasi-free reactions with a spectator proton. The ratio of the number 
of events with a 3rd track to the total number of K^K~ events is 2.9±0.9 % in the mass 
region of 1.50 GeV/c^ < M{nK^) < 1.55 GeV/c^, while it is 8.8±1.2% in the mass region 
of 1.50 GeV/c^ < M{j)K-) < 1.55 GeV/c^, where the A(1520) events dominate. The 3rd 
track ratios in the adjacent regions of 1.55 GeV/c^ < M{nK^) < 1.60 GeV/c^ and 1.55 
GeV/c2 < M{pK-) < 1.60 GeV/c^ are 7.2±1.6 % and 6.4±1.4 %, respectively. A fit to the 
M{nK~^) distribution for events without the 3rd track exclusion shows that the G"^ peak 
height varies only by +0.9 %, while the background level increases by +4.6%. The A(1520) 
peak height decreases by 8.9 % by removing events with a 3rd track. These observations 
indicate the 9"*" peak is likely due to quasi-free •yn reactions. 

In our previous paper, the statistical significance has been estimated from the ratio of 
S/ \/B, where S and B are the numbers of the signal and background events, respectively [l| . 
The toy Monte-Carlo study has shown this method results in large overestimation of the 
significance. The magnitude of the overestimation is still large when using S/ yS + B instead 
of S/^/B. In the current analysis, the S/y/S + B value in the mass region of 1.50 GeV/c^ 
< M{pK^) < 1.55 GeV/c^ is 6.9, which is larger than the significance estimated from the 
A(-21nL) by approximately 2. The ratio of the peak height to its fitting uncertainty gives a 
significance of 5.4 a which is slightly (~5%) higher than that estimated from A(-21nL). 

VIII. CONCLUSIONS AND DISCUSSIONS 

We have observed a narrow peak near 1.53 GeV/ c^ in the nK'^ invariant mass distribution 
from quasi-free 'jn — > K^K~n reactions. The Fermi-motion corrected mass distribution 
is obtained by using the newly developed minimum momentum spectator approximation 
(MMSA). The validity of the MMSA is checked by analyzing the quasi-free 7p i^'^A(1520) 
reactions. The effect of the Fermi-motion on the nK^ invariant mass is studied by using 
the randomized minimum momentum method (RMM), and it has been shown a narrow 
peak with a width much less than ~ 30 MeV/c^ cannot be generated by corrections nor 
selection cuts. The statistical significance of the 6"*" peak has been estimated by a spectrum 
shape analysis using the RMM background functions as well as polynomial functions. The 
statistical significance from the shape analysis is 5.1 a. 
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The G"*" yield is estimated to be 116 ± 21 events in the 1.5 x 10^ K^K~ events. The 
differential cross-section is estimated to be 12 ± 2 nb/sr in the LEPS angular range by 
assuming the isotropic production of the G"^ in the 7n center-of-mass system. 

The Q+/K+K- ratio of (0.8 ± 0.1) x 10"^ is consistent with that of (1.1 ± 0.2) x 10^^ 
obtained by the previous measurement although the detector acceptance is smaller in the 
current experiment due to a longer distance from the target to the spectrometer II]. However, 
the significance of the G"*" contributions in the previous study is highly overestimated because 
it was calculated from the S/ \fB ratio. 

The yield ratio of the G^ to the A(1520) is 0.40±0.09. By considering the partial decay 
branching ratios of r(A(1520) ^ A^Z)/r(A(1520) all) = 0.45 and the acceptance dif- 
ference, the production ratio of the G"^ to the A(1520) is estimated to be 0.15±0.03 in our 
detector acceptance. 

The upper limit on the G"*" production cross-section obtained by the CLAS collaboration 



is 3 nb 



231]. However, due to the different nature of the measurements at CLAS and LEPS, 



the CLAS upper limit is very difficult to compare with the present results, since the re- 
scattering mechanism of the spectator proton (required by CLAS) is unknown and the 
detector acceptances are almost exclusive. Hence we see no conflict between the present 
results and those published by CLAS. 

The LEPS collaboration will analyze new data which were collected with the same de- 
tector setup and an improved luminosity (by a factor of 3). Detailed investigation of the 
angular and energy dependencies of the G"*" photo-production will become possible if the 
peak is confirmed in the new data set. 
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FIG. 1: (a) M{K+K~) distribution for the LD2 runs, (b) MM{-i,K+K-) distribution for the 
LH2 runs, (c) MM{-i,K+K-) distribution for the LD2 runs, (d) MMd{'j, K~)) distribution 
for the LD2 runs. 
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FIG. 2: (a) Pmin distribution for the selected K'^K events, (b) 2-d plot of pp vs. p„ 
non-resonant K^K~ MC events. 
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FIG. 3: MM{'y,K+K-) (left) and MMd{'y,K+K-) (right) distributions for events with 
0.1 GeV/c. 
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FIG. 4: Pmin distributions for events with 2.0 
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FIG. 5: 2-d plots of pmin vs. MM{'y,K ) for simulated signal (G+) events (left) and non-resonant 
K^K~ events (right). 
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FIG. 6: (a) AM = M{nK'^) - MM{j,K-) vs. p^nm with a fit to a 2nd order polynomial 
function (solid curve), (b) M{nK^) (solid histogram) and MM(7, X^) + AM' (dashed histogram) 
distributions for the signal Monte-Carlo events. The dotted histogram is MM{"f^K^). (c) 2-d plot 
of Pram VS. MM{^, K~) for K~^K~ events for LD2 runs. The mean and iblo" are indicated by solid 
and dashed lines, respectively. 



22 





FIG. 7: (a) RMM spectrum for M{nK'^) distribution generated from the signal MC events (solid 
line) and from non-resonant K'^K~ events with 1.50 GeV/c^ < M(ni^+) < 1.55 GeV/c^. (b) 
M{nK^) distribution for (j) events and a fit to the RMM spectra with one seed set (dashed line) 
and three seed sets (solid line), (c) M{nK^) distribution for the sum of non-resonant K^K~ 
MC events and 9"^ MC events and a fit (solid curve) to a mass distribution consisting of RMM 
distributions with three seed regions; (I), (II), and (III). Contributions from each seed region 
is indicated by a dashed histogram, (d) (a) M{nK^) distribution (closed circle) for the MC 
events with a fit to a distribution consisting of RMM spectra and a Gaussian function (solid line). 
The dotted line is the background contribution (the sum of the RMM spectra with fitted weight 
parameters). M{nK'^) distribution for non-resonant K^K~ events (open circle). 
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FIG. 8: (a) 2-d plot of M{K^ K^) vs. Ej-^ and the cut boundary of the (j) exclusion cut (solid 
line), (b) acceptances of the (j) exclusions cuts: the energy dependent cut (closed circle) and the 
constant cut of Mkk > 1-04 GeV/c^ (open circle), (c) M(nK^) distributions for non-resonant 
K^K~ events before the (p exclusion cut (solid histogram) and after the cut (dashed histogram). 
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FIG. 9: The Dalitz plots of M'^{pK') vs. M'^{K+K-) before (left) and after (right) the (j) 
exclusion cut. 




FIG. 10: (a) M{pK~) distribution with a fit to the RMM background spectrum only (dashed line) 
and with a Gaussian function (solid line). The dotted line is the background, (b) MM{^, K^) 
(closed circle) and MM('y,K^) distributions for the LH2 runs. 
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FIG. 11: The Dalitz plots of M'^{nK+) vs. M'^{K+K~) before (left) and after (right) the 
exclusion cut. 
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FIG. 12: (a) M{nK^) distribution with a fit to the RMM background spectrum only (dashed 
line) and with a Gaussian function (solid line). The dotted line is the background, (b) Dalitz plot 
oiM'^{nK+) vs. A'P{pK'). 
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FIG. 13: (a) M{nK+) distribution for events with M{pK-) > 1.55 GeV/c^. A fit to the RMM 
background spectrum only (dashed Hne) and with a Gaussian function (sohd line) in the region 
below 1.65 GeV/c^. The dotted line is the background, (b) The background spectra for the best 
fits to RMM spectra with the wide signal region (dashed line), the narrow region (dotted line), 
and the default region (solid line). 
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FIG. 14: Comparison of the fits with the RMM distributions (soHd Une) and a 2nd-order polyno- 
mial functions (dashed line): (a) in the region of 1.43 GeV/c^ < M{nK^) < 1.65 GeV/c^ without 
the 0+ contribution, (b) with the 0^ contribution, (c) in the region of 1.47 GeV/c^ < M[nK~^) < 
1.65 GeV/c^ without the contribution, (d) with the 6+ contribution. 
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FIG. 15: MM(7, X+) (left) and MM(7, K') (right) distributions with a fit to a mass distribution 
consisting of reversed RMM spectra (sohd line). 
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FIG. 16: M{nK+) (left) and M{pK-) (right) distributions for events with M{K+K-) >1.05 
GeV/c^. The solid lines are fits to the RMM functions plus a Gaussian function. 
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